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We consider exciton effects on current in molecular nanojunctions, using a model comprising a 
two two-level sites bridge connecting free electron reservoirs. Expanding the density operator in the 
many-electron eigenstates of the uncoupled sites, we obtain a 16 x 16 density matrix in the bridge 
subspace whose dynamics is governed by Liuoville equation that takes into account interactions on 
the bridge as well as electron injection and damping to and from the leads. Our consideration can 
be considerably simplified by using the pseudospin description based on the symmetry properties of 
Lie group SU(2). We study the influence of the bias voltage, the Coulomb repulsion and the energy- 
transfer interactions on the steady-state current and in particular focus on the effect of the excitonic 
interaction between bridge sites. Our calculations show that in case of non-interacting electrons 
this interaction leads to reduction in the current at high voltage for a homodimer bridge. In other 
words, we predict the effect of "exciton" blocking. The effect of "exciton" blocking is modified for 
a heterodimer bridge, and disappears for strong Coulomb repulsion at sites. In the latter case the 
exciton type interactions can open new channels for electronic conduction. In particular, in the case 
of strong Coulomb repulsion, conduction exists even when the electronic connectivity does not exist. 

PACS numbers: 73.63.Rt, 73.23. Hk, 73.22.Lp 



I. INTRODUCTION 

Electron transport through molecular wires has been 
under intense theoretical (see e.g. [HQ) an d experimen- 
tal (see e.g. [1,13]) study in the last few years. Theoretical 
studies usually fall into two categories. The first focuses 
on the ab-initio computations of the orbitals relevant for 
the motion of excess char ges through the molecular wire 
while the other floL 1 1 ll] employs generic models to 
gain qualitative understanding of the transport process. 
At the simplest level [H [ll| the wire Hamiltonian is de- 
scribed by a tight-binding chain composed of N sites with 
nearest-neighbor coupling (Huckcl model) that represents 
the electron transfer (tunneling) interactions between ad- 
jacent sites. This model has been generalized to include 
Coulomb interactions between electrons on the same site 
[L2| (Hubbard model) and/or electron-phonon interac- 
tions In the present paper we investigate another 
extension of this model, in which we take into account 
energy transfer interactions between adjacent molecular 
sites. 

Energy-transfer interactions - excitation (deexcitation) 
of a site accompanied by deexcitation (excitation) of an- 
other are well-known in the exciton theory }14l4l6j . In 
particular, Frenkel excitons - neutral excited states in 
which an electron and a hole are placed on the same site 
are readily transferred between sites, and such intersite 
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interactions can accompany the charge transfer processes 
as was shown for charge-transfer excitons [l7| in (quasi-) 
one-dimensional structures [TH, [l| , including polysilanes 
[20l - [22| . The latter show a weak coupling between the 
Frenkel exciton with the admixture of charge transfer 
states and nuclear motions [U . 

In molecular bridges energy-transfer interactions can 
also sometimes have important effects on charge trans- 
fer dynamics. Charge and energy transfer in a lin- 
ear 2,2':6',2"-terpyridine-based trinuclear Ru-(II)-Os(II) 
nanometer-sized array [23j, and one-dimensional en- 
ergy/electron transfer of amylose-encapsulated chain 
chromophores [24j are examples. In addition, it seems 
likely that energy transfer takes place in chemically re- 
sponsive molecular transistors based on a dimer of ter- 
pyridyl molecules combined with ion Co 2+ [25j . 

It should be noted that electron transfer is a tunnel- 
ing process that depends exponentially on the site-site 
distance, while energy transfer is associated with dipolar 
coupling that scales like the inverse cube of this distance, 
and can therefore dominate at larger distances. The im- 
portance of the latter stems also from geometric issues, 
which are related to the dipole-dipole interaction between 
different sites occurring in the vicinity of metal particles 
in molecular nanojunctions. Really, Gersten and Nitzan 
[IE H3 predicted accelerated dipole-dipole energy trans- 
fer near a solid particle (see also 0, [29j]), and in the last 
time a number of works devoted to the exciton-plasmon 
interactions have been published [30l - f33l | that are related 
to physical effects due to the local field enhancement (si 



2 



How will such dipolar interactions affect the conduc- 
tion properties of molecular junctions? This question 
was addressed by Galperin, Nitzan and Ratner by the 
example of a junction composed of one-site- wire and two 
metal leads [40] , where they predicted the existence of 
non-Landauer current induced by the electron-hole exci- 
tations in the leads. To the best of our knowledge, there 
were no analog treatment of simultaneous electron and 
energy transfer (excitons) in multisite bridges. Here we 
address this problem by using the Liouville-von Neumann 
equation (LNE) for the total density operator to derive an 
expression for the conduction of a molecular wire model 
that contains both electron and energy-transfer interac- 
tions. While not a central issue of the present work, we 
note that energy transfer is closely related to heat trans- 
fer through the molecular nanojunction - an issue of im- 
portant consequences for junction stability and integrity. 

Treated separately, the simplest models of exciton and 
electron transport may be represented by tight-binding 
transport models, albeit in different representations. In- 
deed, in the wire Hamiltonian (see Eq.© below), both 
the electron- and energy-transfer terms are binary in 
terms of the annihilation and creation operators for elec- 
trons and excitons, respectively. Their simultaneous 
treatment, however, constitutes a rather complex non- 
linear problem. In this work we combine a tight-binding 
model for electron transport jlft jlu with that of one- 
dimensional Frenkel excitons jl4T - tlo ] to investigate the 
effect of energy transfer interaction on electron transport 
in one-dimensional nanowires. The outline of the pa- 
per is as follows. In Sec|TT]we introduce our model and 
in Sec lIIII we derive a master equation in the eigenbasis 
of many-electron wire Hamiltonian. Sec lIVI is devoted 
to the analytical solution of the problem where we con- 
sider both non-interacting electrons at a site and strong 
Coulomb repulsion at sites. In Sec|V]we show that the 
exciton type interactions can open new channels for elec- 
tronic conduction. In Sec lVII we carry out numerical sim- 
ulations, compare them with the analytical theory and 
show the existence of the "exciton blocking" effect. We 
summarize our results in Sec lVIII In Appendix A we cal- 
culate the eigenbasis of many-electron wire Hamiltonian 
for non-interacting electrons at a site, using the Jordan- 
Wigner transformation [411 ]. In Appendix B we present 
auxiliary calculations. 



II. MODEL 

We consider a spinless model for a molecular wire that 
comprises two interacting sites, each represented by its 
ground, \g), and excited, |e), states positioned between 
two leads represented by free electron reservoirs L and R 
(FigfTJ). The electron reservoirs (leads) are characterized 
by their electronic chemical potentials /i£ and ^.r, where 
the difference ml — /i_R = eVf, s is the imposed voltage bias. 
The corresponding Hamiltonian is 
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FIG. 1: A model for energy-transfer induced effects in molec- 
ular conduction. The right (R = \{r}) ) and left (L — 
manifolds represent two metal leads characterized by electro- 
chemical potentials [ir and hl respectively, each coupled to 
its nearest molecular site. The molecular bridge is a dimer, 
where each site is represented by its ground, \lg) and \2g), 
and excited, |le) and|2e), states. 
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where c^^ (c m /) (m = 1, 2, / = g 1 e) are creation (annihi- 
lation) operators for electrons in the different site states 
of energies e m /, while (c k ) (k G L,R) are creation 
(annihilation) operators for free electrons (energies e k ) 
in the leads L and R. n m f = c^fCmf are the occupa- 
tion operators for the different site states, and site oc- 
cupation operators are given by N m = h mg + n me . The 
operators 6+ = c+ e c mg and b m = c+ lg c me are excitonic 
(creation and annihilation) operators on the molecular 
sites m = 1, 2, while b k > k — c^cy — b^ k , (k, k' e L or R 
) corresponds to electron-hole pairs in the leads. In the 
wire Hamiltonian, Eq. ©, the Af terms represent elec- 
tron hoping between site states of similar energies (i.e. 
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between \g) and between |e) states of adjacent molecu- 
lar sites) , the J terms represent exciton hopping (energy 
transfer) between molecular sites and the U terms corre- 
spond to on-site Coulomb interactions. The molecular- 
leads interactions H conta cts are taken to account for two 
physical processes: V describes electron transfer between 
the molecular bridge and the leads that gives rise to net 
current in the biased junction, while W describes energy 
transfer between the bridge and electron-hole excitations 
in the leads. In JS]) and ([5]) K m is the lead closer the the 
molecular site m {K\ = L, K 2 = R) and H.c. denotes 
Hermitian conjugate. In what follows it will be useful 
also to define the population operators 



A/ = fi2f + n\f 



in the manifolds of ground (/ = g ) and excited ( / = e) 
site levels. 

We consider electronic transport through the molecu- 
lar wire where the leads K = L,R are taken to be each 
in its own equilibrium characterized by its temperature 
T (here taken equal for the two leads) and electronic 
electrochemical potential pk- Therefore, the lead elec- 
trons are described by the equilibrium Fermi functions 
/jf(efc) = [exp((e fc - p, K )/k B T) + l}~ 1 . Consequently ex- 
pectation values for lead operators can be traced back 
to the expression (c^Ck 1 ) = fK{£k)&kk' where 8 k k' is the 
Kronecker delta. The excitonic operators are equal to 



b+ 



The effect of the corresponding interaction 



in the bridge (= HJb^b2 + H.c.) on the charge transport 
properties is the subject of our discussion. 



III. MASTER EQUATION 

Our analysis is based on the LNE, or the generalized 
master equation for the reduced density matrix of the 
molecular subsystem, obtained using a standard proce- 
dure jlo|, Oj], |42| based on taking H conta cts as a per- 
turbation. Briefly, one starts with the LNE for the to- 
tal density operator and use the projectors of the type 
Pxp{t) — picTrKpif) in order to derive an equation 
for the time evolution of the reduced density matrix 
a = TrpTrLp. The calculation is facilitated by invoking 
the so called non-crossing approximation that assumes 
that the effects of different reservoirs (here L, R) and dif- 
ferent relaxation processes (here V, W) are independent 
and additive. This leads to 



da(t) i - 

~¥ TrK I dx[V ' [ym 

where for any operator O, O is the corresponding in- 
teraction representation 



6 mt (-x) = exp[--(H wlre + Hi eads )x] 



OeXp[-(H w ire + Hl ea ds)x] 



(9) 



and where Ttk — Tr^Tr^. 

Consider first terms with the electron transfer interac- 
tions V. Writing the coupling Hamiltonians V n f (Eq.©) 

as 



V„f = c n/ A+ + c+ A nf 



(10) 



( 7 ) where A nf = T,keK n V k {nf) c k , we have V*f(-x) 



S#«(-x)A*$(-x) + c$(-x)A+f(-x) with 

A*>*(-z) = E* g jr B V k {nf) c k eM L h ekx). 

Similarly, writing the coupling Hamiltonian for energy 
transfer W — E„ W n as 



w n = b+e n + b n e; 

(n), 



(11) 



where 9„ = £ fc#fc , 6JfB W U b k'k, then 



W^i-x) = bt int {-x)@^{-x) + b^{-x)&t int {-x) 

(12) 

where 



©r t (-^)= E W^b k>k eMek'-e k )x] (13) 



Bearing in mind that p(t) — a(t)pK where a(t) — 
Tr K p(t) and Eqs.(H0Jl, and ([TZ]). we get for the sec- 
ond term on the RHS of Eq.® 



fi 2 
1 



^Tr K {J dx[V,[V^(-x),p(t)]}} = 



dx{Tr K [VV mt (-x)p K ]a(t) 



V Jo 

- Tr K [V PK a(t)V mt (-x)} - Tr K [V ln \-x)p K a{t)V] 
+ Tr K [p K a(t)V mt (-x)V}} (14) 

In evaluating the RHS of Eq.([TJJ we encounter reservoir 
correlation functions that reflect the reservoir equilibrium 
properties as well as the nature of its interaction with the 
wire. For example, 



Cnf(-x) = Tr K [A n fA+ f (-x)pi 
= E \vt f) \ 2 [l-fK n (e k )]cM-\e k3 

k£K n 



(15) 



Turning to the energy transfer contribution, third term 
on the RHS of Eq.©, we obtain an expression of the 
form ([T4|) with the energy transfer interaction W replac- 
ing V. Using the Wick's theorem, we obtain correlation 
functions of the type 
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D n (-x) = Tr K [®n®Z(-x)p K J 
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w, 



(n) 
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k^k'eK n 



fK n ( £ k)[l ~ fK n { £ k')] expire*: - e k ')x] 

(16) 



Below we get a Markovian master equation in the wide- 
band limit. The full master equation obtained in this way 
constitutes a set of 256 coupled equation for the 16x16 
elements of the wire density matrix, which can be solved 
numerically by diagonalizing the corresponding Liouvil- 
lian matrix. In particular we are interested in the steady 
state solution, ass, which is given by the eigenvector of 
zero eigenvalue. Once ass has been found, the current is 
obtained from 



(7) = Tr(iass) 



(17) 



where the current operator (defined, e.g., as the rate of 
change of electron population on the left of the dashed 
line in Fig. [T]) is given by 



d * ie , - 
e—N = — [H, N] 



* = E 



clc k + hlg + hu 



(18) 



(19) 



keL 



In section [VT] we show some results of such numerical cal- 
culations. To gain better insight of the transport prop- 
erties of this model, analytical simplifications in some 
limits are useful. These are discussed next. 



IV. ANALYTICAL EVALUATION 

It is known [ll[ that for the evaluation of Eqs. © 
and (TT4")) it is essential to work in the representation of 
the eigenstates of the Hamiltonian H W i re + Hi ea( i s that 
defines the zeroth-order time evolution. The use of other 
representations bears the danger of generating artifacts, 
which, for instance, may lead to a violation of fundamen- 
tal equilibrium properties [43j . We thus face the problem 
of diagonalizing a matrix of order 256. This procedure 
may be facilitated by using the pseudospin description 
based on the symmetry properties of Lie group SU(2) as- 
sociated with the two state problem (1/, 2/); / = e,g. 
Such a "donor acceptor" system may be described by the 
"charge transfer" operators = c^ci/ and bf 



c tf c V 



that describe intersite charge transfer 1 — > 2 and 2—^1, 
respectively, in upper and lower states of the molecular 
dimer. The non-diagonal part of H W i re , Eq.([3]), can then 
be written in terms of operators 6/ only 



H 



(nondiag) 



E 

f=9,e 



A f (b++b f )-hJ(b+b g +b+b e ) (20) 



Define also the pseudospin (Bloch) vector in the second 
quantization picture 




bj + bf 
i(b f - b+) 
h 2f - h lf 



f = 9,e 



(21) 



Its components have the following properties: (a) They 
satisfy the same commutation rules as Pauli matrices 
o%2,3 [Hrllll]; (b) the operators A/ = n 2 /+ni/ = 
Em=i,2 c m/ e ">/i / = e,fl (cf. Eq.©) and r{ commute: 
[r{,Xf] = (i = 1,2,3); (c) any linear operator of the 
"donor acceptor" system can be written as linear super- 
position of the operators {r{} and A/. In particular, the 
wire Hamiltonian can be written as 

= ^A e (ei e + e 2e ) + t^ r 3 ( £ 2/ - e lf ) - A f r{] 

f=9,e 

_^( r e r?+r e r | )+ J- U m N m (N m -l) (22) 
m=l,2 

In Eq. (l2"2"T) we have put, without loss of generality, (e\ g + 
e 2ff )/2 = 0. Because the operators A/ and r[ commute, 
A f is conserved under unitary transformations related to 
the diagonalization of H W i re - Therefore, a total 2 4 x 2 4 
space can be partitioned into nine smaller subspaces, i.e. 
the Liouvillian matrix in the required basis is block diag- 
onal with blocks, according to the values of Ay = 0, 1, 2 
(see FigH]): four one-dimensional subspaces for A/ = 0, 2 
for either / = e,g (type I); four two-dimensional sub- 
spaces for A/ = 1 and \f> = 0,2 where / ^ f (type 
II); and one four-dimensional subspace for A e = X g = 1 
(type III). The type I submatrix is diagonal, while four 
state pairs with each pair coupled by the charge trans- 
fer interaction are associated with the four 2x2 blocks 
of the type II subspace. The remaining four states are 
coupled by both the charge transfer and exciton transfer 
interaction and constitute the 4x4 block of subspace 
III. Each of these subspaces is characterized by assigning 
the values (A e , A g ) of total populations in the ground and 
excited states of the two bridge sites. 
Using the identity 



K) 2 = ir{f = (4? = A / -2n 2/ n 1/ = { 

(23) 

the wire Hamiltonian (|22[) can be written in the form 

Hwire = ^A e (ei e + E 2 e) + E U m N m (N m - 1) + 

m=l,2 

+0 For subspaces I 

+ ^3(62/ — Ei/) — A/rf For subspaces II 
+[| E /=ff ,e r f 3 (s2f ~ e lf ) - E /=ff , e A f r{- 
~^f( r i r i + r l r 2)] for subspace III 

(24) 
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FIG. 2: A schematic display of the block structure of the wire 
Hamiltonian. 

This prediagonalization provides an important simpli- 
fication of our problem. From Eqs. (|T7| . (fP5|) . (TH)!) the 
current is given by 

i=j E A /( 6 /- 5 /) = i E A /-2 (25) 

/=S,e /=9,e 

Using Ea.pgJ). this yields 

Z=| E A /^(A/ = 1) (26) 

Obviously A/ = 1 in Eq. (|26[) is another way of saying 
I 



that the current in channel / exists only for the case of 
one of states {/} is occupied and another one of {/} is 
unoccupied. 

Further simplification is made below, when we con- 
sider two specific limiting cases. The first limit, 
U m = 0, describes noninteracting electrons at each 
sites. In the opposite limit with strong on-site Coulomb 
repulsion U m (m = 1,2) is much larger than any 
other energy scale of the problem. In the latter 
case we disregards states with more than one elec- 
tron on any of the molecular site 1 and 2 so we 
need to consider only 9 bridge states: |0i 9 , 2s , 0i e , 2e }, 
|0i g ,0 2ff , li e , he) and |l lfl , l 2g , 0i e , 2e ) in subspaces 

I; |0i g , 2g , 0i e , l 2e ), |0l 9 , Cbg, lie; 02e}) |0l g , llg, Ole, 02 e ) 

and |li g , 2g , 0i e , 2e ) in subspaces II; |li g , 2g , 0i e , l 2e ) 
and |0i g , l 2g , lie, 2e ) in subspace III. 

The diagonalization procedure yields the transforma- 
tion between the eigenstates of the wire Hamiltonian 
and the states of the non-interaction molecular wire, 
\ni g ,n2g,ni e ,n2e) displayed in Figl^J Denoting the col- 
umn vectors of these states by {$} and {x}, respec- 
tively, and the transformation between them by Y", i.e. 
{x} = Y{&}, we can characterized each eigenstate $ by 
the corresponding subspace (A e ,A g ). In this basis, the 
fermionic interaction picture operators (see Eq.©) read, 
for example 



(a\c™}(-x)\f3) = [r + (Ae(a),A s ( a )) x +(Ae(a),A g (a)) c „ /X (Ae(/3),A s (/3))y(Ae(/?),A 5 (/?))]^ 
xexp[i(^(A e (/3) ) A 9 (/3))- J B a (A e (a) ) A 9 (a)))x] 

I 



where (X e (a), X g (a)) denotes the subspace associated 
with the eigenstate a and points to the corresponding 
values of A e and A g , and where (A e (/3), X g (f3j) = (A e (a) + 
l,A 9 (a)) if/ = eand(A e (/3),A g (/3)) = (X e (a),X g (a) + l) 



if f = g. x denotes the transpose matrix x- The relax- 
ation terms in the master equation ([5]) take in this basis 
the forms 



1 f°° 1 

- tf Tr K J dx] y> [V mt {-x),p{t)]U = 2 E T nf{Cnf, aa '<J a 'P'C+ f ^pl2 - f Kn (Ep, - Ep) 

nfa'P' 

- fl< n { E a' - E a )] + C+ faa ,a a ,p,C n f t p,p[f Kn (Ep - Ep, ) + fK n { E a ~ E a ')] ~ {c n f >aa >C+ f a , , f Kn (E a , - Ep,) 

+ c£f,aa'Cn/,a',3'[l ~ fl< n (Ep, - E a ,)}}<Tp,p - <7 aa , {c nf:0l , p, C+ f p, p f Kn {Ep, - E a ,) + 6+ />a , p,Cn f> p, p[l - f Kn {E a - - Ep,)}}} 

(27) 

I 



where 

rn/ = |E Wt f) \ 2 5(e k -e nf ) (28) 

k£K n 
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and 



Tr K dx[W,[W mt {-x),p{t)}] a p 



h 2 ' 

= jE {- B K n [E 0l {\ e + l,X g ) - E a ,(\ e ,\ g + l),/iJf B ][6+ aQ ,6 n ,a'/3'0-/9'/9(*) +ffa/3'(t)6^/ a /6n,a'/3] 
na' [3' 

- B Kn [E a >(\ ei \ g + 1) - Epi(X e + 1, Xg), (J>K n ][b n ,aP'bn f p, a i(Ta' p(t) + &au> (t)b n ,a' P'b*^, p] 

+ bn, aa ^ a 'P'(t)b n ,p'pB Kn [Ep,(X e ,X g + 1) - Ep(X e + l,X g ),(j, Kn ] 
+ b niaa '(T a 'p'(t)b^p,pB Kn [Epi{X e + l,X g ) - Ep(X e ,X g + l),fj, Kn ] 
+ h X, aa > a a' p>{t)b n ,p> pB Kn [E a ,(X e ,X g + 1) - E a (X e + l,Xg),fi Kn } 
P'(tK,p,pBK n [E a ,(X e + l,X g ) - E a (X e , X g + l),n Kn ]} 

I 



(29) 



where 



Bk„ (E a - Ep , [i Kn ) — 



2tt 



k^k'£K„ 



w. 



in) 



kk' 



5{e k - e k > +E a - Ep)f Kn (e k )[l - f Kn (ek')] 



(30) 



In evaluating these forms we have taken the wide band 
limit for the electrodes spectral densities. 

Next consider the diagonalization procedure itself. In 
subspaces I the unitary transformation Y(X e , X g ) is obvi- 
ously the unity matrix. The diagonalization of the block 
matrices in subspaces II and III is carry out in the limit- 
ing cases of zero and infinite on-site interactions. 



A. Zero on site coupling 

The case of zero on site coupling is discussed in Ap- 
pendix A. We find the eigenfunctions and energies of the 
2-site bridge summarized in Table I 



J 



X, = 



Xg=0 



X e = 2 



$(0,0) = 

= |0l s , 2 g, 0l e , 02e) 

E = 

$±(1,0) = 
= ^(|0i ff ,0 2g ,0 le ,l 2e )T 

T|0l fl ,0 2fl ,lle,0 2e )) 

E = e 2e ± A e 

*(1,0)=( jn lff 'o 29 '? le 'H) 
V l U lg; U 2 ff , tie, 2e ; J 

$(2,0) = |0i s ,0 2ff ,li e ,l 2e ) 
E = £ie + e 2e 



Xg = 1 

$+(0,1) = |0i s ,l 2s ,0i e ,0 2e ) 

$_(0,1) = |li s ,0 2g ,0i e ,0 2e ) 
E = 

$(i,i) = y+(i,i)x(i,i) 

E = \{s l e + e 2e )±\J h± 

I | lis 5 ^ff, lie, 2e ) \ 
|llfli 2ff , 0i e , l 2e ) 
|0lg, l 2g , lie, 2e ) 
V |0l s , l 2ff , Ole, l 2e ) / 
$+(2,1) = |0i s ,l 29 ,li e ,l 2e ) 
$-(2,1) = |ll s ,0 2s ,lle,l2e) 



x(M) 



E 



£le + £2e 
I 



Xg = 2 

$(0,2) = |l 9 ,l 2g ,0 e ,0 2e ) 
E = 

$±(1,2) = 

= 7f (|llff; Iff) Ole, le)T 
T|llg, l2 fl , lie, 2e )) 
E = E2e ± Ae 



X(l,2) = 



_ ( |lg, l2g, e , l 2 e) 



|lg, l2g, le, 2e ) 



$(2,2) = |ll 9 ,l 2g ,lle,l 2 e) 
E = Si e + S2e 



where 



/ sin r — cos r — cos r sin r \ 



cos t sin r 
sin r cos t 
cos r — sin r 



sin t cos t 
- cos t — sin r 
sin r — cos r 



(31) 



and where r is given by 

-Jh 

cos 2r = 



^/4A2 + J2ft2 



and sin 2t = 



2A e 



(32) 
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The current in this case is found to be 

(/> = -^A e Im{[a 32 (l, 1) + (74i (1, 1)] cos2r 
n 

+ [a 3 i (1,1)- <7 42 (1,1)] sin 2r- ^ <r_ + (l,A fl )} 

A s =0,2 

(33) 

Indices "+" and "— " in Eq. ([33|) correspond to the func- 
tions $+(l,A g ) and $_(1,A 9 ), respectively, in Table I. 
Indices 1,2,3 and 4 label the the eigenstates of the wire 
Hamiltonian in subspace III. The corresponding ener- 
gies are given by formulas E\ = , E 2 = E--, 
E 3 = E + _, E± = E +t+ where 

E±,± = e e ±±Jh±±y/4A 2 e + J 2 h 2 (34) 



Y + (l, 0) and Y + (0, 1) and the corresponding eigenstates 
and eigenvalues are defined by the same Eqs. (PfBl) . (|T7|) 
and (|4*5|) . respectively, as before (see Appendix A). The 
operator Y + (l, 1) is reduced to (see Appendix B) 

^"-^(-i!) (35) 

Y + (l, 1) is used to obtain the corresponding eigenstates 
<f>(l, 1) = Y + (l, 1)£(1, 1) and eigenvalues £? 1)2 = e e ^Jh. 

Substituting Eq.(|55)l into Eq.JH]) of Appendix A for 
the current, we get for A g = 

(7) = ^A e lm<7_ + (1,0) (36) 
n 

V. CURRENT FROM THE ENERGY 
TRANSFER INTERACTION IN THE WIRE 



B. Rotating-wave approximation 

The calculation of the non-diagonal elements of the 
density matrix <7 a/ g(l, A ff ) in Eg. (f3"3"]) for the current is 
essentially simplified for very weak wire-lead coupling 
when the coherent time-evolution dominates the dynam- 
ics of the wire electrons. This means that the largest 
time-scale of the coherent evolution, given by the smallest 
energy difference, and the dissipative time-scale, deter- 
mined by the electron and energy transfer rates, Y n f and 
Bk„ (E a — Ep , jjLK n ) i respectively, are well separated, i.e. , 
hT nf ,hB Kn (E a - Ep , /J>K n ) < \E a - Ep\ for A e = 1 and 
a ^ (3. Then for A e = 1 and a ^ j3, Eq.© is dominated 
by the first term on the RHS. Consequently, er Q ^(l,A s ) 
can be calculated in the first order of hT n f/(E a — Ep) 
and hB Kn (E a — Ep,fj, Kn )/(E a — Ep). This constitutes 
the essence of a rotating-wave approximation (RWA) [ll[ . 
Within it, one can provide a closed expression for the re- 
duced density matrix elements a a p and for the stationary 
current. We shall use the RWA below in Sec|V] and Ap- 
pendix B. 



C. Strong Coulomb repulsion at sites 

In the limit of strong Coulomb repulsion, U m 
is assumed to be so large that at most one ex- 
cess electron resides on each site. Thus, the avail- 
able Hilbert space for uncoupled sites is reduced to 
three states £(0,0) = |0i ff , 2g , 0u, 2e ), £(2,0) = 
|0i 9 ,0 2fl) li e , he) and £(0,2) = |li 3 , l 2g , 0i e , 2e ) for 

subspaces I; two states £(1,0) = ( l^ 9 ' ® 2g ' ? le ' l 2e \ ) 

\ |Ulg, U 2g , lie, U 2e ) J 

and £(0,1) = ( !? l9 'n 29 'n le '!! 2e ! 1 for subspaces II; 

\ \*-lgi "2g, Ule, U 2e ; J 

and the state £(1,1) = ( \l l9 > ° 29 ' ° le ' l 2e \ ) for the 

\ l 2g , lie, U 2e / / 

now 2-dimensional subspace III. The unitary operators 



In a recent paper [40j Galperin, Nitzan and Rat- 
ner have predicted the existence of non-Landauer cur- 
rent induced by energy transfer interactions between 
a bridge molecule and electron-hole excitations in the 
leads. Here we show that a similar non-Landauer cur- 
rent arises from the exciton type interaction J in the 
wire itself. For simplicity we limit ourselves to elec- 
tron transfer interaction between the wire and the metal 
leads, Eq. (f27|) . and disregard the corresponding excita- 
tion transfer, Eq. (|2l?|) . Also for simplicity we consider 
a large bias limit in the Coulomb blockade case when 
/il > e e and /ifl, < e g , and the states s e ,e g are po- 
sitioned rather far (3> fcsT, H\ J\, |A e |) from the Fermi 
levels of both leads so that fr,(e) = 1 and /r(e) = 
can be taken on the RHS of Eq. ([27|) . Finally, we dis- 
regard electron transfer interaction in the " g" channel, 
i.e. we take A g = 0. Landauer type current would be 
realized in channel "e" when it is isolated from channel 



'<?", i.e. when J = 0, Ti g = T 2g = and A fl 



0. Solv- 



ing Eqs. (|SJ), (|2"T]l in the RWA approximation under these 
conditions and substituting the steady-state solution into 
Eq. (136[) . we get, using also the normalization condition 

EA e =0.1,2 Tr<T ( A ^°) = 1 > 



(I) 



r le r 2 



RWA 



r 



le 



2e 



(37) 



Eq. (l37[) describes the Landauer current and coincides 
with Eq.(21) of Ref . [Ill] (excluding the sign). 

In fact, the current vanishes for Tie = T 2e = even 
when ri ff ,r 2g 7^ 0, since A g = (see Figj3]). Such selec- 
tive coupling to the leads could be obtained for the bridge 
made of a quadruple quantum dot where the lateral ones 
are strongly coupled to the leads [13, |48| . 

Consider now the case when Tie = T 2e = 0; Ti g , T 2g ^ 
0; A g = and J ^ 0. For this case Eqs.©, (HZ]) together 
with Eq.fl3S]) lead to 



(/) = _4eT 9 J a A2 



l-[<7__(0,l) + <r(2 > 0)] 



16A?J2 



H 2 T g 2 J 2 



(38) 
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FIG. 3: A possible physical realization of the selective tun- 
neling configuration, where only \g) levels are coupled to the 
leads. 
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system in state |li 9 , Q^g, 0i e , l2e 
takes place along route such as 02 9 , 0i e , I2 



electron transmission 
1 



|0l 9 , llg, lie, C>2e) 
|0l s , 02 S , 0i e , l2e) 



|0l(,, 02g, lie, 02e) 



llg> 02 ff , Oie, I2 



Step 1 is 

an energy transfer process, steps 2 and 3 rely on T 2g 7^ 
and A e 7^ 0, respectively, and step 4 closes the circle via 
the Tig process. 

Eq. (f3"5]) clearly shows that the current exists only for 
J ^ and A e ^ 0. For small J, (I) ~ J 2 . For large J 
we obtain 



(I) ~ -4er g A; 



l-[(7__(0,l) + a(2,0)] 



16A2 



(39) 



which does not depend on J. In the limit H\J\ , |A e | ^> 
r 9 , Eq.® yields for cj__ (0, 1) = <t(2, 0) = 



(I) 



r 2 „r 



RWA 



2 r 25 + r 



19 



(40) 



In deriving Eq. (|40)) we have not put T\g = T2 g - This 
limit corresponds to the range of validity of the RWA. 
Indeed, it can be shown that Eq. (j40|) can be obtained 
for this model in the RWA (see Appendix C). 

If a (0, 1), cr(2, 0) 7^ 0, the non-Landauer current de- 
creases, since the populations of states |li s , 02 S , 0i e , Cbe) 
{a — (0,1)) and |0i 5 , 2s , lie, he) (<r(2,0)) suppress cur- 
rent. Two latter states are also steady-states in the case 
under consideration (Coulomb blocking, Tie = T2 e = 
A g = 0) along with the states described by FigJU The 
existence of several steady-states corresponds to the pres- 
ence of the respective zero eigenvalues of the relaxation 
matrix. Our numerical calculations give three such zero 
eigenvalues corresponding to three above steady-states. 

If rie,r 2e 7^ 0, state |li g ,0 2g ,0ie,02e) is only steady- 
state that "locks" the current due to Coulomb blocking, 
since A g = 0. Numerical simulations of other situations 
when A g 7^ and non-interacting electrons at a site are 
carried out in the next section. 



VI. NUMERICAL RESULTS 



FIG. 4: Different stages of the energy-transfer induced cur- 
rent, a) energy transfer, er(l, 1) 7^ 0. b) the charge transfer 
to the right lead, c) the intersite charge transfer, er(l,0) 7^ 0; 
the charge transfer from the left lead. 



where for simplicity we put Ti g — T2 g = T g . Eq.(|38|) 
describes a non-Landauer current caused by transport in 
different channels: the intersite transfer occurs in chan- 
nel "e", and the charge transfer between the molecular 
bridge and the leads occurs in channel "g". The inter- 
channel mixing is induced by the energy-transfer term 
J (see FigU]). For example, starting with the molecular 



The results presented in this section are based on di- 
rect numerical solution of Eq. (|8"|). and are in com- 
plete agreement with the analytical solutions when ap- 
plied to the special cases treated in Sections IIVI and [V] 
The numerical solution was carried using the the basis 
of eigenstates of the Hamiltonian H W i re , Eq.Q. Once 
a(t) is obtained from Eq. (JSJ), the expectation value 
of the current is calculated as (/) = Tr(Ia(t)) where 
the current operator was defined by Eq. (|2"oT) . In this 
calculation we have limited ourselves to the case where 
the wire-leads energy transfer coupling W is disregarded 
and, unless otherwise specified, have used the following 
parameters: £i s = e 2g = O.OeV, £i e = £2e — 2.0eV, 
A g = A e = Q.OleV, T lf = T 2f = 0.02eV for / = g, e 
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FIG. 5: The current (I) displayed as a function of the exciton 
coupling parameter J . Vbs = 1.96eV (dashed line); Vbs = 
2.0eV (dotted line); V hs = 4.0eV (solid line) 



FIG. 6: The current (7) as a function of the parameter J in 
the case of non- interacting electrons for Vbs = 4.0eV . The 
current (I) is shown in left panel, and the populations Pi g , 
Pie, Pi a and P2e are shown in the right panel. 



(below we use T to denote the order of magnitude of these 
widths) and T = 100K. The Fermi levels were taken to 
align symmetrically with respect to the energy levels e\ g 
and eie, i.e. [il = (e\ g + e\ e + VLs)/2 and /iji — /il — V& s . 
We also used the value of eA e /h — 2 Ah • 1CP 6 t4 as the 
unit of current (I) . 

Consider first non-interacting electrons. FigsJSJ |5] and 
[7] show the expectation value of the current (J) and one- 
particle populations P n f = Tr{c^fC n f<j) as functions of 
the exciton interaction parameter J. One can see that 
if the imposed voltage bias Vb s is larger than e e — e g , 
the expectation value of the current diminishes when |J| 
increases (we have used J < which is typical to J- 
aggregates, however the trend is similar with J > 0). 
Such a behavior can be understood, using Eg. (|3"4")l for the 
energies in subspaces (III) and Eq. (|33p for the current. 
The latter equation shows two direct contributions to the 
current. The first one has its origin in states of subspace 
(III), the energies of which depend on both A e and J (the 
first and the second terms on the RHS of Eg . (l3"3")) ) . The 
second contribution arises from states of subspaces (II), 
the energies of which depend on A e only (the third terms 
on the RHS of Eq.©). Using Eqs.© with W — and 
(P7)l , the non-diagonal elements of the density matrix on 
the RHS of Eq. (|33f for the steady-state condition can 
be evaluated as o a p ~ e~ 1 -e ( see a ^ so Sec UVB[) . Since 
E- - E+ = -2A e , Ea.(pfe|l. we get for the contribution 
of the third term on the RHS of Eq. fgj} 

yAJm ^-+(l,A ff )~2er (41) 

A g =0,2 

The contribution of the first and the second terms on the 
RHS of Eg. ([55)1 depends on the relation between J and 

A e . 

When h\J\ « A e , Eq.© yields cos2r w 0, sin2r w 
1 , and only the second term in (|33f> gives a contribution to 
the current from the states of subspace (III). Under this 
condition one gets from Eg. (fM)) two doubly-degenerated 
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FIG. 7: Same as Fig[6] for the parameters A 9 = 0, A e = 
O.OleU and V bs = 4.0eV. 

values of energy E\ = E4 = e e + A e and E% — E 3 = e e — 
A e , where the splitting is of the same order of magnitude 
as the hopping matrix element A e . We obtain 

^AJm[a 31 (l, 1) - a 42 (1, 1)] ~ eT (42) 
n 

This contribution is of the same order of magnitude as 
that from the states of subspaces (II). 

In opposite case, h\J\ » A e , cos2r « 1 (again we 
use J < -as in J-aggregates) and sin2r ss 0. In this 
case only the first term in (|33[) contributes to the current 
from the states of subspace (III). For this case we get 
E 2 ~ E4 e e and ~ e e T Jh- This leads to 

^A e Im[ ( x 32 (l, 1) + <7 4 i(l, 1)] - eT^ (43) 

This contribution is much smaller than that of Eqs. (|41j) 
and (|42[) . since the hopping matrix element A e is much 
smaller than the splitting between states 3 and 2, and 
states 4 and 1 due to the exciton interaction (indices 
1,2,3 and 4 label the the eigenstates of the wire Hamil- 
tonian in subspace III) . This can cause the value of the 
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FIG. 8: The current (/) as a function of the Coulomb in- 
teraction parameter Ui — Ui- Vbs = 2.0eV (solid line); 
V bs = 4.0eV (dashed line); V bs = 8eV (dotted+dashed line); 
J = O.OeV (left panel) and J = -0.05eV (right panel). 



FIG. 10: The current (/) plotted against the exciton coupling 
parameter for bias Vbs = 4:.0eV for different energies in the 
e-channel: si e = 1.95eV and e 2e = 2.056^. A 9 = OMeV 
(left panel), A 9 = (right panel). 



l-O.O0eV 
I = -0.02eV 
1= - 0.05eV 



V. [eVl 



J = 0.00 eV 
J =-0.02eV 
J =-0.05eV 



V bs [eVl 



- 0.6 iC 
<" 



V. [eV] 









J=O.O0eV 

■ ■ ■ J = -0.02eV 
- ■ J= - 0.05eV 






J = 0.00eV 

■ ■ ■ J =-0.02eV _ 
J =-0.05eV 










i 




i * 


i 



V h , teV] 



4 « 

< J 



FIG. 9: The current (/) as a function of the bias voltage 
Vbs shown for different values of the exciton coupling pa- 
rameter: J = 0.0 (solid line), J = — Q.02eV (dotted line), 
J = — 0.05eV (dash-dotted line). Left panel - noninteracting 
electrons. Right panel U\,U2 = oo. 



total current to decrease. In other words, the transitions 
3 — s- 2 and 4 —> 1 do not participate in electron transfer 
due to their large splitting for h\J\ » A e . This is in a 
sense " exciton blocking" of electron transmission through 
the bridge. 

Next we turn to situations where electron-electron in- 
teraction is taken into account. Fig|S] shows the current 
(I) as a function of the Coulomb interaction parameter 
U\ = U2- Fig|9] depicts the current (I) as a function 
of the bias voltage Vb s for different values of the exci- 
ton coupling J for the case of non-interacting electrons 
as well as for the case of infinite on-site interaction be- 
tween electrons. The "exciton blocking" effect seen for 
non- interacting electrons (smaller current for larger \J\) 
disappears in the case of Coulomb blocking. 

This is supported by Eq. (f3l>| that does not show 
a direct contribution of the states of subspace (III) 
to the current, and the above evaluation of the term 
^A e Imer |_ (1, 0) . The point is that in the case of in- 
teracting electrons, subspace (III) includes only states, 



FIG. 11: Same as Fig. [9] except that ei e = 1.95eV and e% e = 
2.05eV\ 



which are acted upon exciton interaction (see Figj2]). 
Moreover, in the case of Coulomb blocking, the effect 
of exciton- induced current exists (Sec fV|) . 

The effect of "exciton blocking" depends on the en- 
ergy detuning £2/ — £\j in channel " /" for a heterodimcr 
bridge. Figs HO] and [11] show the current (/) as a function 
of J for £2 e — £ie = O.leV. (I) is seen to increase for small 
|J|, then to decrease as \ J\ becomes larger. This can be 
related to the modification of resonance conditions when 

Finally, figures [121 [151 [Til and [151 show more of the 
system behavior for the model with U\,U2 — 00. FigfT2l 
shows the current as a function of |«7| for different values 
of the imposed voltage bias Vbs- If Vbs is large com- 
pared to the energy difference between the excited and 
ground site energies, the current behaves in accordance 
with Eq. (f3"51) . If Vbs is close to this energy difference, the 
current increases initially with \J\, and then decreases to 
zero. Furthermore, in accordance with Eq. ([55|) . the left 
panel of FigfT3lshows that the steady-state current is zero 

for the initial condition a (0, 1) = 1. The steady-state 

current is zero also for the initial condition <r(0, 0) = 1, 
since the latter state relaxes to a (0, 1) = 1. Figures [TH 



11 




-0.06 -0.04 

J [eVJ 



^ 11 










P l5 




1 , 1 


... 




\ 




\ 


- - P le 




~~ \ 






\ 






•s. _ „ ■ ■*" 






/ 

/ 









Time [fs] 



FIG. 12: Current as a function of |J| for the initial population 
of state |0i 9 , 2g , lie, 2e > equal to 1. A 9 = 0.0, A e = O.OleV, 
ri s = r 2g = 0.02eV, r ie = r 2e = O. V bs = 2.0eV (solid 
line), Vbs = 4.0eV (circles - numerical simulations, dashed - 
calculations with Eg. (1381) ). 




Time [fs] Time [fs] 

FIG. 13: Current (I) as a function of time for different initially 
populated many-electron states: |0i 9 , 2ff , 0i e , 2e ) - solid, 
|li 9 ,02 9 ,0i e ,0 2 e) - dashed, J = -0.05eV and V bs = 8.0eV. 
Left panel: A 9 = 0, A e = O.OleV, Fi 9 = T 29 = 0.02eV, T le = 
F 2e = 0, |0i 9 , 29 , 0ie, lae) - dot-dashed, |0i 9 , la 9 , 0i e , 2e ) 
- squares, |0i 9 , 2g , lie, 2e ) - dotted. Right panel: A 9 = 
O.OleV, A e = 0, Fi 9 = T 29 = 0, T le = T 2e = 0.02eV, 
|0i 9 ,0a 9 , li e ,0 2e ) - dot-dashed, |0i 9 , l 2g , 0i e , 2e > - dotted, 
|0i 9 ,0 29 ,0i e , l 2e ) - squares. 

and [15] show the time dependence of the current and one- 
particle populations for different initial conditions corre- 
sponding to the absence of relaxation in e-channel and 
^-channel, respectively. 

VII. CONCLUSION 

We have developed a theory of electron transport 
through a molecular wire in the presence of the effect 
of dipolar energy-transfer interaction between the sites 
in the wire. We found that such interaction, which leads 
to exciton excitations in the wire, cannot in general be 
disregarded. We used a model comprising a two two- 
level sites bridge connecting free electron reservoirs. Ex- 



FIG. 14: Current and one-particle populations P n f = 
Tr{c^fCnf>a) as functions of time for the initial popula- 
tion of state |0i 9 , 2g , 0ie, l 2e ) equal to 1. J = — 0.05eV, 
V bs = 8.0eV, A 9 = 0.0, A e = O.OleV, Ti 9 = T 29 = O.OleV, 

rjv/.ie = r 2e = o.o. 
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FIG. 15: Current and one-particle populations P n j = 
Tr(c^fC n f/cr) as functions of time for the initial popula- 
tion of state |0i 9 , 29 , li e , 2e ) equal to 1. J = — 0.05eV, 
Vbs = 8.0eV, A 9 = 0.0, A e = O.OleV, F lg = T 29 = 0, 
Fie = r 2e = 0.02eV. 

panding the density operator in the many-electron eigen- 
states of the uncoupled sites, we obtain a 16 x 16 density 
matrix in the bridge subspace whose dynamics is gov- 
erned by Liuoville equation that takes into account in- 
teractions on the bridge as well as electron injection and 
damping to and from the leads. Our consideration has 
been considerably simplified by using the pseudospin de- 
scription based on the symmetry properties of Lie group 
SU(2). We studied the influence of the bias voltage, the 
Coulomb repulsion and the energy-transfer interactions 
on the steady-state current and in particular focus on the 
effect of the excitonic interaction between bridge sites. 
Our calculations show that in the case of non-interacting 
electrons this interaction leads to reduction in the cur- 
rent at high voltage for a homodimer bridge. This ef- 
fect can be called "exciton" blocking. The effect of "ex- 
citon" blocking is modified for a heterodimer bridge, and 
disappears for strong Coulomb repulsion at sites. In the 
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latter case the exciton type interactions can open new 
channels for electronic conduction. In particular, in the 
case of strong Coulomb repulsion, conduction exists even 
when the electronic connectivity as defined above does 
not exist. 

To end this discussion we note that in this work we 
have investigated a molecular bridge connecting metal 
leads. It is worthy to note that the geometry considered 
could modify the effect of dipolar e nerg y-transfer inter- 
action between the sites in the wire |26|. This issue will 
be considered elsewhere. 
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VIII. APPENDIX A. NON-INTERACTING 
ELECTRONS AT A SITE 



The unitary transformations Y(X e ,X g ) = I for sub- 
spaces (I). As to subspaces (II), Hamiltonian correspond- 
ing to the second line of the RHS of Eq. (IM1) where A/ = 1 

Xf (/' 7^ /) can be diagonalized, using the unitary 
transformation 





(r[ 


(3) 


= f f \ r{ 






where 






cos2z9/ = 




sm2tf/ = 



cos2tf/ -sin 2??/ 

1 
sin 2-d f cos2i?/ 



S 2 f - Slf 



(£2/-£i/) 2 +4Aj 
( £2/ - £l/ )2 + 4A2 



(45) 



The matrix elements of T* are connected with the uni- 
tary transformations Y(X e , X g ) for subspaces (II) by for- 
mula = (l/2)Tr(& n Y+ajY) where a n and tjj are 
Pauli matrices. 



_ / cos i9f sintf/ 



defined by 



Y+(X f = l;X f , = 0,2) - , 

which enables us to obtain eigenstates 

$ + (A/ = l;V=0,2)\ - + 
$_(A/ = 1;A// =0,2) )- Y W- 1 ^/' 

xx(A / = l;A / / =0,2) 
and eigenvalues 

E ± (X f =l;X f > =0,2) = i[A e (e le +e 2e ) + (e 2 /-ei/) 



(46) 

0,2) 
(47) 



± \ Ke 2 f -£i/) 2 



-4A 2 ] 



(48) 



for subspaces (II). Here the many-electron eigen- 
states of the uncoupled sites are given by x(l,0) = 
/ |0 lg ,0 2ff ,0 le ,l 2e ) \ , . _ ( |0 lg ,l 23 ,0 le ,0 2e ) \ 

^ |Oi s ,o 2ff ,i le ,o 2e ) J' X[ ^ L> { |ii s ,o 2s ,o le ,o 2e ) y 

x(i, 2 ) = ( and ^ - 

\ I J-lff) i-2g, i-le, V2e) J 
( |0lg, l 2 g, lie, l2e) j 
\ 2 g, lie, l2e) / 

Taking the expectation value of the current, Eq. ([2"ol . 
we get 



(I) 



2c: 



{ A / Im<7_ + (A / =l;A / 0-X] Im<7 ^( 1 ' 1 ) 



X fl =0,2;/ 



cu/3 




x [Y+(l, l)x+(l, l)(A B 6 e + A g b g )x(l, l)Y(l, \)U) 

(49) 

where we put — i? 2 for Xf = 1 and Xf — 0,2 that 
follows from Eq.(P| and used (ifc£(A/ = 1; Xf = 0, 2)) = 
Tr{a 2 (r) = 2Imcr_+(A / = l;A/< = 0,2). Indices "+" and 
" — " in Eq. p9|) correspond to the functions $+(1, X g ) and 
<J>_(1,A S ), respectively, in Table I. 



B. Unitary transformation for subspace (III) 

The calculation of Y + (l, 1) is more involved. Consider 
for brevity a homodimer bridge with e ng — 0, e ne = 
e e and A g = 0. Bearing in mind future generalizations 
of our model to A-sites, we shall transform the Paulion 
operators (&t ,bf) to fermion operators (fit, j3f) through 
the Jordan- Wigner transformation [4ll |43, [50J : 



A. Unitary transformations for subspaces (II) 

Consider subspaces (II). In the limit U m — 0, the 
matrix T Eq. (|4"4")) . with matrix elements T£j = 

(l/2)Tr[a n Y+(X f = 1; Xf = 0,2)a 3 Y(X f = l;X f , = 
0,2)] describes a rotation by mixing angle 2§f around 

axis "y". Y(Xf = l;Xf = 0,2) is an unitary operator 



Pe = b e, Pt = bt,!3g = cxp (iirb + b e ) b g) 
0+ = b+ exp (-iirb+be) (50) 

Then H W i re , Eqs.Q and (|20p . can be rewritten for sub- 
space (III) in terms of the fermion operators as 

H wi re(K = X g = l)= e e ~A e (/3 + +/3 e )-hJ(f3 + p g +p+p e ) 

(51) 
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Eg. (l5lT) is a quadratic in Fermi operators and can be 
diagonalized in two stages. Its "excitonic" part H ex = 
—hJ((3ePg + P^fle) is readily transformed to satisfy the 
condition H ex = . HejOjdj if we take (49j 

a j = y^iPg sin ^ + /3 e sin ^), 

ej = -2Jcos^, j = 1,2 (52) 

where a-, are also Fermi operators. The correspond- 
ing occupation number basis set contains 2 2 = 4 eigen- 
functions of the system. The si ngle- excited states are 

given by o+|0) = E^/l/> = VWls) + (-l)''" 1 ^) 
/ 

where |0) = |li g , 2ff , li e , 2e ) is the "vacuum" state, 
and \g) = \0 lg , l 2g , lie, 2e ) and |e) = |l l9 , 2g , 0i e , l 2e ) 
are the states with the corresponding donor acceptor pair 
excited. The eigenstate with two excitations can be writ- 
ten down in terms of the Slater determinant 



Substituting Eq. (l3"TT) into Eq. (|49| for the current, we 
get Eq.® for A g = 0. 



IX. APPENDIX B. UNITARY 
TRANSFORMATION FOR SUBSPACE (III) FOR 
INTERACTING ELECTRONS AT A SITE 

In the limit of strong Coulomb repulsion, the oper- 
ator F+(l, 1) is reduced to that defined by Eg. (1331) in 
accordance with Egs. ([521) . since the "hopping" operator 
Fj = (— l)- J (A e /v / 2)(a^" + a,j) has no nonzero matrix el- 
ements involving states with a single excitation |e) and 
\g) (see Appendix A). 

Substituting Eq. (l35l) into Eq. (|49|) for the current, we 
get Eq.([Ml) for A g =0. 



44J°) 



t>j 2 e 



"319 



\eg) = -l(-iy> - (-iyi]\eg) 

(53) 



with energy ei + e 2 = equal to that of the vacuum state 
where \eg) = |0i g , l 2ff , 0i e , l 2e ). The wire Hamiltonian 

A„ = 



can be written down in terms of aj as H W i re (X, 



j 3 ' J 



Fi = 



I) = e e + J2j Hj where Hj = F 3 + (-l)^Ja 

(— l) J (A e /v / 2)(aj h + aj) is the "hopping" operator with 
the only nonzero matrix elements involving states which 
differ by a single excitation: (0|.F>+|0) = (-l) j A e /V2, 



\a j2 F n a+a+ 



10} 



(— I) 31 A e /V2. The eigenstates and 
eigenvalues of H W i re (X e = \ g = 1) can be calculated now 
as follows. #(1,1) = r+(l,l)x(l,l) where F+(l, 1) is 

/ |0) 



given by Eqs.([SI|) and (JSU), £(1, 1) = 




$(1,1) 



/(|0) + |e 5 ))sinr-(| e ) 
J_ (| e ) + |. g ))sinr + (|0)- 
y/2 (Ie> -1.9)) cost + (|0) - 

\(|0)-|e fl ))coBT-(|e) 

l*i) 

1*2) 
1*3) 
1*4) 



, and 



f- |y))cosT 
I eg)) cos t 
\eg))smT 

~ l5 , ))sinr 



(54) 



X. APPENDIX C 

The steady-state solution of Eqs. ©, ([27]) in the RWA 
approximation gives for the case under consideration 



ct(0,0) = a(0, 2) = <7 ++ (0, 1) = a+^(0, 1) = a_+(0, 1) = 

(55) 

and (7 (0,1) and a(2,Q) are arbitrary. Putting 

cr__(0, 1) =cr(2,0) =0, we get 



<7_+(l,0) 



-8A, 



■{r M ,i s rm(l,0) + T M>2B Tra(l, 1)} 

(56) 

and Tm(l,0) = (T Mt2g /T M ,Xg)Tra(l, 1). Then using 
the normalization condition 



Tm(l, 0) + Tm(l, 1) = 1 



and Eq.®, we obtain Eq.gO]) of SecEl 



(57) 
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